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Abstract 

The functional flow equations for the Legendre effective action, with respect to 
changes in a smooth cutoff, are approximated by a derivative expansion; no other 
approximation is made. This results in a set of coupled non-linear differential 
equations. The corresponding differential equations for a fixed point action have 
at most a countable number of solutions that are well defined for all values of the 
field. We apply the technique to the fixed points of one-component real scalar field 
theory in three dimensions. Only two non-singular solutions are found: the gaussian 
fixed point and an approximation to the Wilson fixed point. The latter is used to 
compute critical exponents, by carrying the approximation to second order. The 
results appear to converge rapidly. 
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Our aim is to find a reliable and accurate analytic approximation method of general 
applicability to non-perturbative quantum field theory [2]-[4]. We believe that such 
a method must be based, in practice, on the Exact Renormalisation Group 1 [1] 
(ERG), for reasons explained in ref.[2]. 2 We particularly consider here the non- 
perturbative low energy behaviour - the Wilson fixed point [1] - of three dimensional 
one-component (Z 2 invariant) scalar field theory, which is also the continuum limit 
of the three dimensional Ising model and physically the universality class of critical 
binary fluids and liquid- vapour transitions in classical fluids [11]. We choose to study 
this because of its relative simplicity but the method - just a derivative expansion of 
the Legendre effective action combined with the ERG - is certainly, appropriately 
modified[10], of much wider applicability to non-perturbative quantum field theory. 
In this letter we concentrate on developing a practical framework for such an 
expansion. The explicit calculations are only taken to second order, which is not 
sufficient (with our choice of cutoff) for our estimates of exponents to be competitive 
with the worlds best present estimates [11]; nevertheless, already the differences, 
between our results and those, are 6 or less times their quoted error, and our error 
appears to roughly halve with each new order of approximation. Finally we should 
emphasize that, in calculating directly the effective lagrangian, we may derive much 
more of physical interest than just critical exponents. 

The reader may well be asking why such a natural approximation scheme, built 
around ideas nearing twenty years old[l], has not been well studied before. It 
is as well to realise however that there are many contenders for methods of 
approximation, that sound perhaps equally natural, but do not work. It is also 
true that the computation must be organised carefully: merely expanding vertices 
in invariant polynomial combinations of momenta will soon grind to a halt through 
shear complexity. We list below methods that do not work in general. 

a) A sequence of truncations of Dyson- Schwinger equations, and other methods where 
renormalisability is a serious problem [2]. 

b) A sequence of truncations of the ERG, because these do not converge beyond a 
certain point, and because no completely reliable method could be found to reject 

1 a.k.a. Wilson's renormalisation group 

2 For examples of other recent work on the ERG see refs.[5]-[10]. 
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the many spurious solutions that are also generated [3] -[5] (see however [6]). 

c) Direct numerical approaches to (truncations of) the ERG [2]. 

d) Approximations other than those that can be formulated for one particle irreducible 
parts of the vertices of the Wilson effective action [2]. 

e) Momentum (i.e. derivative) expansion of the ERG with smooth cutoff, if the 
effective width over which the cutoff varies is any less than the cutoff itself [2] . 

f) Momentum scale expansion of the ERG with sharp cutoff because only certain 
truncations are calculable in closed form[4] - but this is not good enough, c.f. (b). 

Because of (d) we use flow equations formulated directly for the Legendre effective 
action with respect to changes in some infrared cutoff scale A. This gives directly 
the one particle irreducible parts of the Wilson effective action as detailed in ref.[2]. 3 
Because of (c) we must approximate the momentum dependence and the only 
efficacious method we know of, is a momentum expansion[2]. Because of (e) and 
(f) we use a smooth cutoff with an effective width at least A (that is with little 
variation on scales less than A). This means that the cutoff will only somewhat 
suppress, and not eliminate, low energy modes. Finally because of (b) we perform 
only the derivative expansion and do not at all expand in the field <p. 

The partition function is defined as 



The notation is the same as previously [2] -[4], so two-point functions are often 
regarded as matrices in position or momentum (q) space, one-point functions as 
vectors, and contractions indicated by a dot. We work in D euclidean dimensions 
with a single real scalar field (p. The definition differs only in that here we mean 
Sa to be the full bare action for the theory; the C = C(q,A) term is taken to 
be an 'additive cutoff', for convenience. Previously we took a formulation with 

3 It should be realized that all the ERG equations are straightforwardly equivalent 
to each other, up to choice of cutoff function; for example the Wegner- Houghton 
equations[12] are the Polchinski equations[13] in the limit of sharp cutoff, while these 
are equivalent to Wilson's equations [1] under the momentum dependent change of 
variables (p — > (py/~K that eliminates the cutoff function K [13] from the propagator. 
For other comparisons see refs. [2] and [14]. 




(1) 
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a multiplicative cutoff 6> e , so that a cutoff massless propagator took the form 
9 e {q, A)/q 2 . Of course the two are trivially related, in this case 



£ (q, A) = q 2 C(q, A)/[l + q 2 C(q, A)] 



(2) 



The reason for the change is that the smooth cutoff will have to scale anomalously 
along with ip to reach non-gaussian fixed points, and this scaling is homogeneous for 
smooth additive cutoffs, but takes the form (2) for smooth multiplicative cutoffs. 

For C(q, A) to be an infrared cutoff it follows that C should be small for q < A, ideally 
C tending to zero as q — > 0, and q 2 C(q, A) should be large for q > A, tending to 
infinity as q — > oo. The ultra-violet regularisation need not be discussed explicitly [2], 
in fact since we will concentrate exclusively on the neighbourhood of fixed points it 
can be ignored entirely. From (1) we have 

5W dC' 1 5W 



+ tr 



DC- 1 5 2 W \ 
' 8.16.1 ) 



^8J dA 8 J \ dA 
which on rewriting in terms of the Legendre effective action V gives (as in ref.[2]), 



1 

-2 tr 



\_d_c_ ( gr_V 

C dA'{ + '5ip5ipl 



T is defined by T[(p] + \ip.C~ x .<p = —W[J] + J.(p, where now ip = 5W/5J is the 
classical field. As in ref.[2] it is helpful to write the trace as an integral over 
momentum space and factor out the .D-dimensional solid angle: 



dA = -2h dq 



g^- 1 dC(g,A) 
dA 



1 + C. 



S 2 T 



(q, -q) 



(3) 



C(q, A) 

where O = 2/[T(D /2)(4tt) d / 2 ] is the solid angle of a (D - l)-sphere divided by 
(27r) D , the brackets (• • •) mean an average over all directions of the momentum q, 
and we have written the momentum dependence of the traced inverse two-point 
function explicitly. 

At a fixed point the field scales anomalously as <p ~ A d ^ with d v = \{D -2 + 7/), 
r\ being the anomalous scaling dimension. Dimensional analysis then indicates that 
we must have 

C(q, A) — > A v ~ 2 C(q 2 / A 2 ) (4) 
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for some C as we approach the fixed point, if V is to become independent of A 
- i.e. fixed - in that limit. What is happening is that the anomalous scaling of 
(p induces also an anomalous scaling of the width over which the smooth cutoff 
changes; it is because this width is zero in the sharp case that anomalous scaling 
may be incorporated in that case automatically [12] [4]. Now we rewrite the equation 
in terms of dimensionless quantities 4 via (4), q i— > Aq, <^(Aq) h- > A d ^~ D (p(q) : and 



t = ln(A /A). Finally we rescale the fields and effective action as (p t— > ip^Q/2(, 
r i— > (Q/2QF to absorb the factor 0/2 in (3). £ is a normalization factor which will 
be chosen later for convenience. The result is 



(- + d ip A ip + A d -D)T[ l p} = 



C dqg - 1 
Jo 



9 »c(f) +2 _; 



C(q 2 ) dq 



(5) 



(q, 



In here A v = ip.-^ is the 'phi-ness' counting operator: it counts the number of 
occurences of the field <p in a given vertex. Ag may be expressed as 



^ f dp . . 
D + J(2^ IP{P)P 



(2n) D dp» 5(p(p) 

i.e. the momentum scale counting operator +D. Operating on a given vertex it 
counts the total number of derivatives acting on the fields <p. 

We write V as the space-time integral of an effective lagrangian expanded in 
powers of derivatives. Each linearly independent (under integration by parts) scalar 
combination of differentiated fields <p will be accompanied by its own general (t 
dependent) coefficient function of the fields: 

r>] = J d D x {v(<p, t) + ±(d^) 2 K(^ t) + {d^fH^ t) + (n v ) 2 H 2 ( v , *) + •••} . 

(6) 

(We have not written all 3 four-derivative terms). Substituting this into (5), 
expanding the Right Hand Side (RHS) of (5) up to some maximum power of 
derivatives - this is the approximation - and equating coefficients of the linearly 
independent basis set, yields a set of coupled non-linear differential equations for 
the coefficient functions. 



4 From now on we drop the tilde on the scaled C. 
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The derivative expansion of the RHS of (5) may be computed as follows. Regard 
1 



1 + C. 



1 + C. 



as a differential operator: 

S 2 T r 1 

6ipd<p> ^' 



q)= ld D xd D ye- 1 ^ 
d D x Q 



1 + C. 



S 2 T 

dipdip 



(x,y) 



where 



= p -^qx 



Q = e 



1 + C. 



S 2 T 

dipdip 



-i -l 



,zq.x 



(7) 
(8) 



Q is a function of q, and v?( x ) an d its derivatives evaluated at x. Using (8), we see 
that Q satisfies the differential equation 



+ 



l+V l+V 



VQ 



-zq.x 



C(-n) 



S 2 T 

dipdip 



3 zq.x 



Q 



(9) 



where C and £ J are also regarded as differential operators - acting on all terms 
to their right. For convenience the function V(</?(x), q) is introduced and defined to 
be the expression obtained by dropping from 

5 2 T 



C(q 2 ) e~ zc l- x 



dipdip 



.zq.x 



(10) 



all terms containing differentials of (p. With this definition the derivative expansion 
of Q may be straightforwardly performed by iteration of (9). To complete the 
computation of the RHS of (5), it only remains to perform the angular average 
and compute the resulting radial q integrals. But the angular averages are easily 
translated into invariant tensors: (q^q u ) = q 2 b^ v j D etc. 

We will require that t-independent (fixed point) solutions for the coefficient func- 
tions: V((p), K((p), Hi(ip), H 2 ((p), ... be non-singular for all real (p. In addition we 
will require V to be parametrized so that K(0) ^ 0. 5 As a result it will be seen that 
the differential equations have at most a countable number of fixed point solutions 
(at least for the first two orders and presumably to any order of approximation 
in the derivative expansion). The coefficient functions of linearized perturbations 
about these solutions - ie. of the eigenoperator spectrum - will be required to grow 



The obvious stability requirements do not need to be separately addressed since 
unstable solutions are not allowed by the equations. 



no faster than a power of ip as ip — > oo. It will be seen that this forces power law 
growth (as <p — > oo), the power being given by the scaling dimension, and results in 
a quantized spectrum. 

For a general cutoff function C(q 2 ), for example based on an exponential, this 
would appear to leave rj as an undetermined parameter however. In this case 
there is presumably a range of reasonable approximations, which can be thought 
of as labelled by the effective width of the cutoff (e.g. via (2) modified by the 
value of K(0)). On the other hand, dimensional analysis indicates that if C(q 2 ) is 
homogeneous in q 2 then the flow equation (5) is invariant under a global scaling 
symmetry. This symmetry then overdetermines the equations so that solutions exist 
only for certain discrete values of 77. (This is similar to the case of sharp cutoff [4] [12] 
where the flow equations are invariant under the phi-ness scaling symmetry <p — > \ip). 
We pursue this form for the cutoff here, and therefore require C(q 2 ) = q 2k for k a 
non-negative integer. (The proportionality constant may be set to 1 by a rescaling of 
(p) . The scaling symmetry is characterised by the following dimensions, as trivially 
follows from (6) and the definition of V (see above (3)): 

P/J = [?/*] = 1 > [<p] = k + D/2, 

(11) 

hence [V]=D, [K] = —2(k + 1), [H l } = -D - A{k + 1) . 

(Of course these dimensions should in no way be confused with their physical scaling 
dimensions). 

Using (2) one finds [1 - 6 e (q, A)]/q 2 = l/[q 2 (l + q 2k+2 )]. This is a massless 
propagator with effective U.V. cutoff as it appears for the Wilson action[2]. One sees 
that this choice of cutoff corresponds to a form of Pauli-Villars or higher derivative 
regularisation. It may be that other cutoffs give stronger convergence, so the more 
general case deserves its own investigation. From our general arguments [2] we 
expect that the momentum expansion has slower convergence the higher the value 
of k. On the other hand we require k > D/2 — 1 if the momentum integral in (5) 
is to be U.V. convergent to all orders in the momentum expansion. Therefore we 
choose k to be the least integer larger than D/2 — 1. The cutoff function is now 
completely determined. The inequality k > D/2 — 1 may be derived by thinking 
graphically about the RHS of (5): At a given order p 2m (m > 1) in the momentum 
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expansion one obtains an expansion in one loop graphs with N = 1, 2, • • • vertices 
that behave as ~ q 2m as q — > oo, AT + 1 propagators that behave as ~ l/q 2m as 
q — > oo, and an insertion ~1/C(q). Putting these together with the measure factor 
~ g D one finds the stated inequality. (Actually one should consider also the m = 
case for which propagators go as ~ 1/q 2 and vertices are constant, but this yields a 
weaker constraint.) 

We are finally ready to consider the explicit example. Thus from now on we set 
D = 3 and k = 1. To lowest order we drop all the derivative terms (d^ip etc.) from 
the RHS of (5). This means that the coefficient functions K, Hi , H2 , . . . satisfy 
linear equations given by the vanishing of the LHS of (5). This implies that at 
fixed points K((p) oc ip~ 2r, / 1+v , while the requirement that K(ip) be non-singular 
and K(0) 7^ implies we must have 77 = (as expected at this order) and K{ip) 
a constant, which we set to 1 by (11). Similarly we determine Hi(ip), • • • = 0. Q 
is just 1/(1 + V) to this order, and by (10), V = q 2 [q 2 + V"(<p)]. (We use primes 
to denote differentiation w.r.t. cp). Substituting (7) into (5) and performing the q 
integral we find 

J^(V?, t) + \<pV'(tp, t) - 3V(<p, t) = -\j y/2 + V"(<p, t) . (12) 

Here, and from now on, we fix Q = l/2n. Solutions to this equation are qualitatively 
very similar to those of the equivalent sharp cutoff equation[7] [3]. For any putative 
non-trivial fixed point solution such that V(ip, t) = V(ip), some analysis shows that 
either it exists for all real <p, in which case it behaves, for <p — > 00, as 

^^ v + 37iW + 0(lM) (13) 

for some positive constant A (an isolated one parameter set, as follows from the 
analysis of (15)), or else it is only defined for \cp\ < <p c because the solution ends at 
a singularity of the form 

/ 9 \2/3 

V((p) ~ \J^) (^c - V?) 4/3 for ip w +v? c . 

(We have not shown the lower order singular behaviour and regular part. Note 
that neither here nor at 0(d 2 ) does V(<p) diverge at the singularity.) There is 

7 



also an analogous "slow rollover" behaviour [3]). With dV/dt = 0, (12) is a second 
order ordinary differential equation. Thus at a fixed point, solutions are generically 
labelled by two parameters; however there are two conditions to be met: that V((p) 
have the behaviour (13) and that V((p) be symmetric under ip <-> —<p, equivalently 
V'(0) = 0. ( In a non-symmetric theory this condition is replaced by requiring 
(13) also for (p — > — oo with a possibly different constant A). Thus there are at 
most a discrete set of acceptable fixed point solutions to (12). We have searched 
numerically for these and find only two: the trivial gaussian solution V(cp) = l/y/18 
and the solution "<9°" displayed in fig.l, which we expect to be an approximation 
to the Wilson fixed point. It may be characterised by (12) and the "effective mass" 
a = V"(0) = —.5346, negative as expected[3]. 6 For small perturbations about this 
solution we write V(<p, t) = V(<p) + SV(<p, £), with, by separation of variables, 

6V{<p,t) oc v(ip)e Xt . (14) 

Substituting in (12) one obtains 

yv\ V ) + {\-Z)v{ V ) = \v"{y)[2 + V"{<p)r m • (15) 

Again we expect solutions labelled by two parameters, however in this case by 
linearity we can choose v(0) = 1, and by symmetry we require i>'(0) = 0. Thus 
solutions v(<p) are unique, given A. However for large <p, v(<p) is generically a 
superposition of ~ </? 6-2A , which is that expected from (14) and scaling arguments, 
and ~ exp(i[30A] 3 / 2 </? 8 ). Requiring zero coefficient for the latter restricts the 
allowed values of A to a discrete set. We found just one positive eigenvalue, which 
yields the correlation length critical exponent [1] through v = 1/A, and determined 
the exponent of the first correction to scaling uj = — A, where A is the least negative 
eigenvalue. The corresponding solutions v((p) are displayed in figs. 3. Our results 
are displayed in Table 1. These can be compared to the equivalent sharp cutoff [7] [3] 
results v = .6895, uj = .5952. (It is worth remarking that with these methods one 
very readily improves on the results for exponents given in these refs.) 

All numerical values reported in this letter have been determined to an accuracy 
greater than the number of significant figures displayed. 
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Working to 0(d 2 ) in the derivative expansion we determine, as before, that 
coefficient functions Hi ((f) : ■ ■ ■ vanish identically. To lowest order Q is 1/(1 + V) 
with V = q 2 [q 2 K((p) + V"(<p)}, but now we must iterate (9) to 0{d 2 ) using C(-U) = 
-□ and ^ = V" - K'd^d^ - \(d^) 2 K" - (Oip)K' - KU. After a long but 
straightforward computation, involving also integration by parts of (7) and the 
angular average, and performing the q integrals, the RHS of (5) may be cast in the 
form (6). Comparing both sides of (5) we obtain 

W + hl + ^V-3V=- _ 1 ,-" /4 _ (16) 



dt 



2 y/Ky/v" + 2y/K 



dK I,, x t„ „ (, V\ f 1 24KK" -19(K') 2 

and — + -(l + r))<pK' +nK = 1--W ; V- ' , (17) 

Ot 2 V 1 V A \48 K S 2 (V" + 2^Kf 2 



1 58V" K'VK + 57 \K') 2 + (V"') 2 K 5 (V"') 2 K + 2V"'K'VK + (K 1 ) 



l\2 



+ 



48 K(V" + 2VK) 5 / 2 12 VK(V" + 2VK) 7 / 2 j 

Considering again the fixed point equations (with dV/dt = 0, dK/dt = 0) we 
see that, by substitution of (16) into (17), it is possible to cast these equations 
in the form of two simultaneous second order ordinary differential equations. A 
little thought shows that this generalises: at any order of a derivative expansion, 
substitution of the lower order equations in the higher order equations will reduce 
them to simultaneous second order differential equations, arising ultimately as a 
consequence of the term £ £ in (5). There are now two ways the solutions can 
end at singularities, 7 but we find the essential point remains the same: if V(ip) and 
K((p) exist for all real cp, and are not constant, then again they must behave, for 
(p — > 00, to leading order according to their scaling dimension: 

VUf) = Ay <// 1+ " + ; 1+V -2(1-V)/1+V + O ( -(6~2V)/1+V) (18) 

4,/6(5 - v)A K A v V J 



K{<p) = A K ^n+v + I (1 71/2)2 p-V+nVi+ri + o U- 8 /^) (19) 

6(1 + 7^/6(5 -?7)A y v ' 



for some positive constants A V ,A K . (Now an isolated two parameter set; this 
then accounts for the four behaviours expected. We have also investigated (16) 
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V(ip) ~ const. (ip c — </?) m , K(ip) ~ const. (ip c — (p) n with (m, n) = (4/3,2/3) and 
(8/5,-4/5). 
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and (17) numerically, and have checked analytically that (18) and (19) are the only 
possible power law behaviours (as (p — > oo) excepting several special possibilities 
which however overdetermine the five parameters). 

Since rj is so far a free parameter, solutions are generically labelled by five 
parameters. However there are now also five conditions to be met - which may 
be given as (18), (19), V'(0) = 0, K'(0) = 0, and K(0) = 1. This last condition can 
be chosen using the scaling symmetry (11). Thus for these "non-linear eigenvalue 
equations" there are at most a discrete set of allowed values for rj and associated 
solutions V(cp), K(<f). We have searched numerically for these and again find only 
two: the gaussian solution (rj = 0, K{<p) = 1, V((p) = l/\/l8), and an approximation 
to the Wilson fixed point displayed in figs.l and 2. It is characterised by (16), (17), 
the solution rj given in Table 1, and a = V"(0) = -.3782. 

As before, we linearize about this fixed point solution: 8V(cp,t) = ev(cp)e^, 
SK((p, t) = e k(cp) e - ^, e infinitessimal. Linearity and symmetry allow us to choose 
v(0) = 1, v'(0) = and k'(0) = 0, while requiring that v and k behave for large tp 
according to their scaling dimension as <£> 2 ( 3-A )/ 1 + ? ? and <£> -2 ( A +'?)/ 1 + ? 7 respectively, 
and not as growing (and oscillating 8 ) exponentials of powers of <p, once again turns 
these equations into eigenvalue equations for A. As before we find just one positive 
eigenvalue and determine the first negative eigenvalue. These yield the values for v 
and lo given in Table 1 below. The corresponding solutions for v and k are given in 
figs. 3 and 4. 



Approx' 


V 


V 


UJ 


o(d°) 
o(d 2 ) 

Worlds Best 




.05393 
.035(3) 


.6604 
.6181 
.631(2) 


.6285 
.8975 
.80(4) 



Table 1. Results from the first two orders of the derivative expansion compared 
to the worlds best determinations [11]. (The latter are combined results from the e 
expansion, summed perturbation theory and lattice methods). 9 All other exponents 
pertinent to the fixed point itself are related to rj and v by scaling relations. 

8 ~ exp{c(l ± iy/3)ip p }, where p = | and c > is a calculable constant. 

9 A recent lattice Monte Carlo RG study however gives a much lower rj, and u and v 
much closer to our 0(d 2 ) results [15]. 
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However, we also find solutions v and k with eigenvalue A = 0. This is a simple 
consequence of the scaling symmetry (11). These are not physical but rather are 
components of a redundant operator [14] corresponding to a reparametrization of 
the effective action: In general a change of variables ip x = (p x + e $> x [0\ in (1) (with 
J. if replaced by J.(p) induces a change in the effective action of 5T = F.j£ with 
F x [(p] =e exp(—W[J])$ x [5/8J]exp(W[J]) and a change to the definition of the 
cutoff functional [14], in effect for the classical field 5ip x = F x . A general choice of F 
that leaves (6) invariant in form at this order is F x [f] oc {/(</?(x)) +ai fl 9 ( j^(x)} for 
/ any function and a any constant. A redundant operator corresponding to this F 
must satisfy (v, k) oc (fV — 3aV , fK' + 2f'K — aK) for some choice of / and a. 
Such is the case for the solutions at A = 0, with /(<£>) = 5ct</?/2 and a^O. 
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Fig.3. The v(ip) component of the relevant scaling operator (associated with u), 
and the least irrelevant scaling operator (associated with a;), at order d° and d 2 . 
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Fig.4. The k(cp) component of the relevant scaling operator (associated with z/), 
multiplied by a factor of 10 (for display purposes), and the least irrelevant scaling 
operator (associated with u) : at order d 2 . 
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